Understanding the structure of the first atomic contact in Gold 
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We have studied experimentally the phenomena of jump-to-contact (JC) and jump-out-of-contact 

(JOC) in gold electrodes. JC can be observed at the first contact when the two metals approach 

each other while JOC occurs in the last contact before breaking. When the indentation depth 

rT) between the electrodes is limited to a certain value of conductance, a highly reproducible behaviour 

y—{ in the evolution of the conductance can be obtained for hundreds of cycles of formation and rupture. 

f^ Molecular dynamics simulations of this process show how the two metallic electrodes are shaped into 

fSJ tips of a well-defined crystallographic structure formed through a mechanical annealing mechanism. 

. We report a detailed analysis of the atomic configurations obtained before contact and rupture 

fS I of these stable structures and obtained their conductance using first-principles quantum transport 

.^fH calculations. These results help us understand the values of conductance obtained experimentally 

^ in the JC and JOC phenomena and improve our understanding of atomic-sized contacts and the 

t evolution of their structural characteristics. 
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Background 



Metallic atomic sized contacts can be created by Scanning Tunneling Microscopy (STM) [H [2] or by mechanically 

controlled break junctions (MCBJ) [Tl|3]. In such nanocontacts the electrical conductance is closely related to their 

minimum cross section. Therefore, by recording the conductance while the electrodes are displaced with respect to 

Q each other (traces of conductance), one can infer the atomic structure of these contacts. However, to understand the 

O structures formed at the contact it is necessary to make use of theoretical models. Uzi Landman et al. |3] pioneered 

the use of molecular dynamics (MD) simulations to follow the variation of the minimum cross section during the 

^ process of stretching of a nanocontact. Later, Untiedt et al. j5], by experimentally studying the phenomena of jump- 

J> to-contact (JC) in gold and combining MD and electronic transport calculations, were able to identify the formation 

]^ of three basic structures before contact between the two electrodes, although a limited analysis on the conductance 

values was presented there. 

Trouwborst et al. [5] have also studied the phenomena of JC and JOC using indentation loops where the maximum 

conductance was limited to IGq, where Go = 2^ (quantum of conductance). These experiments showed that the 

^^ elasticity of the two electrodes is one of the relevant parameters to explain these phenomena. Despite the above, there 

(T^ is not yet a unique picture that correlates the experiments with the MD and transport calculations, regarding the 

^H different atomic structures that can be found at the contact. 

•^ On the other hand, experiments, together with molecular dynamics and electronic transport calculations based on 

• '-j density functional theory, show how very stable structures can be obtained by repeated indentation. This has been 

rS described as a mechanical annealing phenomenon ^T^. Limiting the maximum conductance value (5Go for gold) in the 

^ process of formation and rupture of a nanocontact leads to reproducible and atomically sharp pyramidal electrodes. 

This technique has recently been used by other authors [8] to prepare tips in situ for low temperature STM. 

In this paper we show experimental results of the JC and JOC phenomena for gold, analyzed simultaneously. 
We study the most probable configurations before the formation and breaking of nanocontacts with pyramidal form 
obtained from MD simulations emulating the process of mechanical annealing. As found earlier [Sj, the contacts 
can be classified into monomer, dimer and double contact. In order to correlate with the experimentally obtained 
conductance values, we calculated the conductance of these structures using first-principles quantum transport models. 



Methods 

We have used a STM, where the tip and sample were two gold electrodes with purity 99.999%. The experiments 
were done at 4.2K and cryogenic vacuum atmosphere. In order to obtain the conductance of the contacts the electrical 
current was measured while applying 100 mV constant bias voltage between the gold structures. Figure lA shows 
traces of conductance in a gold nanocontact, measured in units of Gq during the process of formation (red) and 
rupture (green). Insets show some snapshots from our molecular dynamics simulations. These correspond to the 
initial structure (top figure) and the final structures before breaking (bottom right) and just after contact formation 
(bottom left). Figure IB is a zoomed area around IGq of figure lA, where the phenomena of JC and JOC can be 
clearly observed. In order to quantify the jump occurring in these two processes we define two conductance values for 
JC {Ga, Gb) and two values for JOC {Gc Gd)- These values correspond to the conductance values before and after the 
jump. We have performed thousands of indentations and recorded the values of these points. Representing Gb vs Ga 
for the JC case and Gd vs Gc for JOC, we can obtain a colour density plot as shown in Figure IC for JC and in Figure 
ID for JOC. Lighter colours are less probable values than darker colours. To emulate the movement of the STM 
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FIG. 1: How to build a density plot. On the top left hand side we show a typical trace of conductance of gold at 4K during 
the process of breaking (red) and forming (green) a contact. The top right hand size is a zoom near IGo to define the values 
before and after the jump-to-contact, Ga and Gb, and jump-out-of-contact, Gc and Gd- The bottom figures show colour density 
plots where dark colours represent those values of conductance that appear more frequently (left for JC and right for JOC). 



and simulate the tip and surface annealed mechanically, we used MD simulations with embedded atom potentials. 
Density function theory (DFT) based calculations are performed to obtain the electronic transport in the simulated 
structures [9]. For the MD simulations we have selected an embedded atom potential [10] because elasticity of the 
electrodes seems to be one of the key parameters in the processes to be studied [6] , and these empirical potentials are 
fitted to reproduce the experimental elastic properties of bulk materials. 

Furthermore the computational cost with this simulation method is low, which makes it an appealing tool since we 
need to simulate tens of these cycles of breaking and formation of the nanocontact. 



Using MD we have analyzed the same structures described in detail in reference [7], but now we focused on the 
type of contact formed. The two initial configurations of the nanocontacts are shown in Figure 2. Structure A is 
built with 525 gold atoms. This initial structure is stretched until the contact is broken by displacing the two top 
and bottom atom layers (represented in blue in the figure). After breaking, the direction of the displacement of these 
layers is reversed so that the two sides are brought together until contact. The temperature in the simulations is 
4:.2K. In this case the temperature is scaled in every cycle of breaking and formation of the contact. The indentation 
process continues until the minimum cross-section formed has 15 atoms. Then, the whole cycle starts again, breaking 
and forming the contact for a total of 20 cycles (see moviel of supplementary material on ref. |7j). The second 
structure studied (B) is shown in Figure 2, composed of 2804 gold atoms. In this case the indentation is limited to 
cross-sections of 25 and 15 atoms (movie2 and movieS at supplementary material on ref. |7]). The temperature here 
is kept constant and equal to 4:.2K during the whole simulation by scaling the velocities of all atoms every time step 
(every femtosecond). The strain rates applied are between 10® and 10^°s~^, typical of MD simulations [TT]. Note that 
the ratio of length of the contact to minimum cross-section is very different in these two structures (5 for structure A 
and 2 for structure B), therefore exploring a system with a long and narrow constriction and another of a short and 
wide nanowire. As shown previously [7] , structure A reaches a stable configuration formed by two pyramidal tips after 
repeated indentations. This configuration is formed after cycle 11 and it remains stable for the following 9 cycles. In 
each of these cycles, although the pyramidal shape remains, there are differences in the atomic configurations right at 
the contact, as shown in Figure 3. These are the configurations we study and describe in this paper in detail. For the 
case of structure B, because of the initial shape, the formation of the two pyramidal tips occurs from the very first 
cycle and, again, only differences are observed in the very last atomic configuration forming the contact. We have 
performed electronic transport calculations based on DFT [HI H^ for both structures, A and B. These calculations 
have been carried out with the help of our code ANT.G, which is part of the package ALACANT [Q' and implements 
the non-equilibrium Green's function formalism as a module to the popular code GAUSSIAN09 [12J. Here, due to 
the large number of atoms, we have employed a very basic basis set consisting only of one 6s orbital and one electron 
the remaining 78 electrons being part of the pseudopotential. 




structure A 



FIG. 2: Structure and configurations of contacts. The two initial configurations used in the MD simulations are shown; 
structure A, long and narrow contact and structure B, short and wide contact. 



Results and discussion 

Experimental results of the JC and JOC in gold are shown in Figure 1 and Table I. Figures IC and D show the color 
density plots obtained for gold when representing Gb vs Ga for the case of JC, and Gd and Gc for the case of JOC. 
Note the presence of two very distinct areas in the JC plot, corresponding to configurations with a high probability. 
In the case of JOC we can distinguish clearly one area of high probability. More details about these experiments are 
presented in reference [5] . For clarity, we include in Table I those pairs of conductance that appear more frequently 
in the experimental measurements. We should mention that for all traces studied in gold the phenomena of JC or 
JOC are always observed, unlike in other metals[5]. For JC we observe three pairs of values that occur with higher 
frequency which we named as maximum 1, 2 and 3. In the JC maxima 1 and 2 correspond to jumping from a value 
of O.OlGo to a value of 0.9400 and from a value of 0.05Go to 0.986*0. These two peaks are easily observed in Figure 
IC as one large area of high probability. The last maximum corresponds to a jump from 0.09Go to 1.77Go, which is 
the second spot shown in Figure IC. On the other hand, on breaking the nanocontact, only two maxima have been 
identified, one where the contact breaks for conductance values of 0.92Go, clearly seen in Figure ID, and another one 
when it breaks at conductance values of I.6OG0, which appears very faint in the figure. Note that these two values 
are close to those obtained for the first and third maximum in the JC case. 



TABLE I: Experimental Values of conductance that appear more frequently in the case JC and JOC. Pairs of values {Ga,Gb) 
and {Gc,Gd) for JC and JC-JOC respectively, that appear more frequently in the density plot of Figure 1. 



Pairs of values obtained in the density plots in fig.l 


Phenomena 


1 maximum 

{Ga, GhjGo 


2 maximum 

{Ga, Gb)Go 


3 maximum 

{Ga, Gb)Go 


JC 


(0.03,0.94) 


(0.05,0.98) 


(0.09,1.77) 


JOC 


(0.01,0.92) 


- 


(0.01, 1.60) 



As mentioned above, we make use of molecular dynamics simulations and DFT calculations of conductance to 
understand these experimental measurements and observations. Figure 2 shows the two structures studied using MD, 
as described earlier. In Figure 3 we show some snapshots of the configurations found just after the contact between 
the two tips and just before breaking a nanocontact. Three basic atomic structures are found: a monomer and a 
dimer (Figures 3a and 3b), and a double contact (Figures 3c-3e). For the case of a double contact we have identified 
different geometries, three of which are shown in this figure. We introduce, for first time, the concept of a double 
dimeric (Figures 3c and 3d) and monomeric (Figure 3e) contact. We define a double dimeric contact as the one where 
the contact is between two atoms facing two other atoms, while we define a doble monomeric contact as a contact 
where two atoms are contacting each other. Another interesting point is that for the double dimeric contact we have 
identified two possible structures, one where two atoms are perpendicular to the other two (Figure 3c), which we call 
transversal configuration (D.C. Dimeric T), and one where two atoms are parallel to the other two (Figure 3d), which 
we call parallel configuration (D.C. Dimeric P). 
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FIG. 3: Structures are the point of contact or before breaking from MD simulations. Representative configurations obtained 
from MD simulations right before contact or right before breaking are shown: (a) dimer (b) monomer (c) double contact dimeric 
transversal (d) double contact dimeric parallel and (e) double contact monomeric. 



Table II shows the probability of finding a monomer, a dimer, or a double contact (all possible configurations 
for D.C) in the MD simulations right before contact and right after contact, for the 2 initial structures and different 
indentations. Note the limited statistics in these results, since only 10 cycles have been computed for the first structure 
and 9 cycles for the second one. Nevertheless, we can see some interesting results. 

For the case of structure A, with a large ratio of length to minimum cross-section, we observe that the most probable 
configuration both at JC and at JOC is a dimer. The monomer and the double contact have similar probabilities. 
This result is in agreement with reference |13| . The situation for the second structure, B, with a small ratio of length 
to minimum cross-section, is significantly different. In this case, when the indentation between the two tips is limited 
to 15 atoms in cross-section, the configuration at the contact is the same in all cycles, a double contact, although 
we observe the formation of the different double contacts described in Figure 3c-e. Clearly, very stable pyramidal 
structures are formed in this case. The robustness of the tip imposes the repetition of a certain kind of structures. 
When the indentation between the two tips increases to a value of 25 atoms in cross-section we should note that the 
traces do not repeat between cycles, and therefore, different structures are formed. In this case, for JC, the double 



TABLE II: MD results of the type of first contact or last contact (JC or JOC) in the structures A and B annealed mechanically. 
This table contains the probabilities of finding a dimer, monomer or double contact in the MD simulations for the two structures 
studied A and B during the formation or the breaking of the contact and for different indentation values (15 atoms or 25 atoms 
in the minimum cross-section). 



Percentage of cases of type Monomer, Dimer and D.C. 


Phenomena 


Monomer 


Dimer 


Double Contact 


JC Structure A 15 iden 


20 


60 


20 


JOC Structure A 15 iden 


30 


60 


10 


JC Structure B 15 iden 








100 


JOC Structure B 15 iden 








100 


JC Structure B 25 iden 


22 





78 


JOC Structure B 25 iden 


22 


56 


22 



contact is still predominant, while for the JOC the probabilities have the same trend as in structure A (dimer being 
the most probable). 

In order to correlate the results from molecular dynamics to the experimental measurements it is necessary to 
calculate the conductance of these atomic structures. 

TABLE III: Electronic conductance calculated by DFT on typical contacts obtained from MD structures. 



Structure and Value of Conductance Go 


Metal 


Dimer 


Monomer 


Doble Contact 


Au 


0.92 ± 0.07 


0.97 ± 0.15 


1.73 ± 0.02 



Table III shows the values of condctance obtained from electronic transport calculations based on DFT for the typical 
first or last contacts proposed: monomer, dimer and doble contacts. The table includes the values of conductance 
obtained with their standard deviation. We can observe that the monomer values of conductance are in the range 
1.20 — 0.76Go, with an average value of 0.97Go. That is because during the process of rupture and formation the 
monomer can be localized closer or further away from the rest of the contact. Another important factor that can 
change the conductance of a monomer is the total number of neighboring atoms to the central atom in the contact, 
which can be different while remaining a monomer structure. Both factors are responsible for the spread in the 
conductance values of a monomer. On the other hand, the deviations in the conductance values for dimer or double 
contact structures are significantly smaller, around 0.07Go and 0.02Go respectively, the average conductance value 
being 0.92Go for a dimer and 1.73Go for a double contact. These results indicate that, on average, dimer and 
monomers have similar values of conductance while double contacts have significantly larger conductance values. It 
seems clear then that the maximum obtained experimentally for JC and JOC, with conductance values of 1.77Go and 
I.6G0 respectively, (maximum 3 for JC and maximum 2 for JOC in Table I) correspond to the formation of a double 
contact. The results for the other maxima obtained experimentally are not so clear since the average conductance 
values obtained for a monomer and a dimer in the calculations are very similar. This seems to indicate that the 
two first maxima obtained experimentally in the JC must correspond to configurations in a dimer and in a monomer 
geometry. According to MD simulations, the most likely configuration both in JC and JOC is a dimer (except in 
special cases of very stable tips) , although monomers can also be formed. 

Figure 4 combines the experimental results with the calculations of conductance in order to identify the experimental 
values of JC and JOC with structures detected by MD and calculated by DFT. The center column between the 
experimental density plots of JC and JOC indicates the average value of conductance obtained from the simulations for 
each geometry (double contact, monomer and dimer). The thickness of the rectangles around each geometry indicates 
the standard deviation. It is clear from this plot that the top high frequency events in the density plots corresponds to 
a double contact and the bottom high frequency events corresponds to monomer and dimer configurations. Although, 
as we mentioned above, it is difficult to distinguish the monomer and dimer using our theoretical model, we can see 
that the average of conductance of monomers is above the one of the dimers. If we add to this that we would expect 
a higher tunnel conductance (on average) prior to the formation of a monomer we can label maximum 1 and 2 as 
dimer and monomer respectively. 
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FIG. 4: JC and JOG density plots together with conductance calculations of different geometries of the contact. Inside the 
experimental density plots we have marked the average conductance values after or before the jump as obtained from DFT 
electronic transport calculations with their deviations. 



Conclusions 



Experiments of JC- JOC show that certain structures are more Hkely to occur than others. This depends on the metal 
and on the process of breaking/formation and the type of structure at the electrodes. Simulations and calculations (MD 
and DFT) of these experiments show that three basic atomic structures are formed at the contact: monomers, dimcrs 
and double contacts. We have identified within the double contact structure several different atomic arrangements 
that we named double dimeric contact (parallel and perpendicular), and doble monomeric contact. According to DFT 
electronic transport calculations, double contacts have an average value of conductance of 1.73Go, which correlate 
very well with one of the peaks observed experimentally both for JC and for JOC. This configuration is also obtained 
in JC and JOC from the MD simulations and, for some very stable tips, is the dominant configuration. Monomers 
and dimers, however, are difficult to distinguish from the simulations since their average conductance values are very 
similar (0.97Go and 0.926*0 respectively). In the case of JOC these two peaks can not be resolved. Interestingly, 
the conductance values are somehow lower than in the case of JC, which could indicate the most likely formation of 
stretched contacts. 
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